Understanding the Heavy Fermion Phenomenology from Microscopic Model 
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We solve the 3D periodic Anderson model using a two impurity cluster DMFT. We obtain the 
temperature v.s. hybridization phase diagram. Approaching the quantum critical point (QCP) both 
the Neel and lattice Kondo temperatures decrease and they do not cross at the lowest temperature 
we reached. While strong ferromagnetic spin fluctuation on the Kondo side is observed, our result 
suggests the critical static spin susceptibility is local in space at the QCP. We observe in the crossover 
region logarithmic temperature dependence in the specific heat coefficient and spin susceptibility. 
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Heavy Fermion phenomenon is among the most in- 
tensively studied subjects in condensed matter physics 
. Experimental information accumulated over the past 
thirty years has revealed many unconventional aspects of 
the heavy Fermion physics [llQl- 

Heavy Fermion physics is derived from a local moment 
band hybridizing with an extended conduction band. 
The hybridization induces the competing Kondo (Jk) 
and RKKY interactions 0. At Jk ^ T, the physics 
is dominated by those of the two bands separately with 
the hybridization as a perturbation. As the temperature 
is lowered, depending on the strength of Jk, different 
physics may develop. In the region where Jk ^ W (the 
conduction bandwidth), the RKKY interaction prevails. 
The RKKY is a long range exchange mediated by the 
conduction electrons near the Fermi surface and oscillates 
with k = 2kF asymptotically At T < J^/W, the RKKY 
interaction can induce a transition to a magnetically or- 
dered phase. In the crossover regime on this side, the 
Kondo behavior, though subdominant, would still show 
up in various measurables due to its non-analyticity in 
terms of the energy cutoff, e.g. temperature. As Jk 
is increased, the Kondo effect becomes more important 
and eventually dominates. Here the Kondo screening be- 
gins at high temperatures T ^ Jk where the conduction 
electrons near the Fermi energy starts to screen the lo- 
cal f-moments. If this remains so as the temperature is 
lowered, there would not be enough conduction electrons 
to completely screen the f- moment lattice Actually 
in heavy Fermions the entire conduction Fermi sea gets 
involved in screening. As the temperature is lowered, 
conduction electrons farther away from the Fermi surface 
participate in the Kondo screening. The system may then 
be described, to the leading order, as a band of local mo- 
ments whose magnitude is progressively reduced. These 
reduced local moments still hybridize with the conduc- 
tion electrons which have not participated in the Kondo 
screening and live near the conduction band bottom. 
Hence, on the Kondo side, the RKKY correlation, as T 
is lowered, becomes more ferromagnetic (FM). The FM 
spin correlation remains at further lower temperatures 
when the heavy Fermi liquid is formed. Macroscopically, 
the FM behavior is related to the lattice Kondo energy, 



To, which is proportional to the Fermi energy. Tq can 
be defined in terms of the saturated homogeneous spin 
susceptibility in the Fermi liquid phase, Xg^o ~ ^/^Oi 
where C is the Curie constant. Should Tq approach zero, 
strong FM spin fluctuation would be observed. 

It turns out that the thermodynamics related to the 
continuous condensation of the local moments into the 
heavy fermion fluid is quite universal |15|, as contrasted 
with that more material specific in the low temperature 
region. At low temperatures, various phases may de- 
velop, including a superconducting phase I']. One in- 
teresting possibility is that the competing Kondo and 
RKKY interactions result in a quantum phase tran- 
sition (QPT) without the interference of any other 
phases. Such a situation is observed experimentally in 
CeCui^,:Au^ HQ and YbRh2Si2 0. 

Two different scenarios have been proposed for the 
heavy Fermion QPT. One is the Hertz-Millis-Moriya the- 
ory 9] , which is applicable when the energy scale of the 
Kondo screening is much higher than that of the mag- 
netic ordering near the quantum critical point (QCP). 
As a result, the local moments are fully screened in the 
quantum critical regime. However, by comparing the pre- 
dicted critical exponents with the experiments, the the- 
ory is found to be unapplicable in many cases jli 0. ■ 
In the second scenario, one expects the magnetic order- 
ing (T/v) and lattice Kondo screening (Tq) energies vanish 
simultaneously at the QCP 0i (^^e also Hflj). Since 
the critical spin fluctuations at fc = 2kF and fc = are 
associated with the T/v and vanishing Tq, respectively, 
it is interesting to know what kind of critical magnetic 
mode may develop in the neighborhood of the QCP. We 
will show later that near the QCP the two critical modes 
strongly interact with each other and the critical spin 
fluctuation becomes local in space (see also [2]). In this 
scenario, the local moments survive at the QCP. 

It is important that the heavy Fermion physics as we 
understand be obtained from a microscopic model. To 
this end the main theoretical difficulty lies in treating on 
equal footing the Kondo and RKKY interactions. Many 
bosonic mean field theories, like the Hertz-Millis-Moriya 
and slave-boson theories, fail because they rely 
on the order parameter of either the magnetic or Kondo 
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phase and miss the properties of the other. On the other 
hand, a fermionic mean field theory, Hke the dynamical 
mean field theory (DMFT) , would allow the possible 
orders to develop and compete and is more desirable. 

The DMFT extends the Weiss mean field theory to de- 
scribe fermions. The single impurity DMFT was applied 
to both the Kondo and magnetic ally ordered phases in 
heavy Fermions (see also Ref.0]). Besides the over 
estimation of the Neel temperature due to the lack of 
the magnetic fluctuation, this approach can not capture 
properly the renormalization of the RKKY interaction. 
A partial solution to this problem is to add an RKKY in- 
teraction to the model and allow the renormalization of 
the RKKY by extending the Weiss approximation to the 
interaction 15]. The resulted formalism, the so-called 
extended DMFT (EDMFT), is able to describe quahta- 
tively the heavy Fermions 16] in both the Kondo and 
antiferromagnetic (AFM) phases. However, it predicts a 
first order phase transition due to the local mean field 
treatment of the RKKY interaction 0, 0| . 

A parallel path in studying the heavy Fermions is 
through the two impurity problem. This model contains 
the dynamics of both the Kondo and RKKY interactions 
and is solvable 0,^3- '^^^ shown, at the particle- 
hole symmetry, there was a non-Fermi liquid fixed point 
separating the Kondo and magnetic phases. However, it 
is difficult to extend the properties of a multi-impurity 
model to a lattice of impurities. 

In this letter, we combine the two impurity model with 
DMFT so that a lattice of impurities can be described. 
This overcomes the difficulties of both the single impurity 
DMFT and multi-impurity approaches. It handles the 
Kondo and RKKY interactions in a more balanced way. 

We consider the periodic Anderson model in 3D 
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(1) 

with — —^^i^i cos ki. We divide the lattice into 
two interpenetrating sublattices, A and B. The unit cell 
is then doubled. Applying the cavity method we 
obtain an effective local action: 



The Weiss function Qq^ is determined self-consistently 
as follows. First, we use quantum Monte Carlo method 
(QMC) and obtain the impurity Green's function. 
Then, from the Dyson equation for the impurities, we 
get the impurity self-energy [E™^']xy- The lattice self- 
energy is constructed to be: 
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from which we obtain the local Green's function 
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with /(fc) = exp(ikx)- By identifying the local (within 
a unit cell) Green's function on the lattice with that 
of the impurity model, we form a self-consistent loop 
[13. l2ll| . While solving the impurity model we can 

dcf 

measure the z-direction spin susceptibility xxyIt) — 
{Tr SiiX,T)Si{Y,0)), with SiiX) - n^^^ - n{^^. The 
lattice susceptibility is obtained in the same way as the 
lattice self-energy given in Eq.©. 
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FIG. 1: The calculated phase diagram. The two lines, 
Tn{V) and To{V), do not cross at the lowest temperatures 
we reached. The inset shows the low temperature saturation 
of the static homogeneous spin susceptibility a,t V = 0.28. 
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We study the phase diagram of temperature v.s. V 
at fixed U — 1.2 and Ef ~ —0.15. To avoid crossing 
the band gap, we change the chemical potential /i along 
with V so that the free {U — 0) particle density per 
site at T = is fixed at Nf^l"" = 2.5423. The result- 
ing physical density changes slightly as V increases and 
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is always greater than and close to the half- filling |22l |. 
We study the AFM and paramagetic (PM) phases and 
the transition between them. In solving the two impurity 
problem we use the QMC. We always use U At ^ 1 where 
At — f3/L and L is the number of time slices in QMC. In 
each DMFT iteration, we perform QMC sweeps ^ 10^. 
Away from the phase transitions around 10 DMFT iter- 
ations are usually enough to converge the results. Near 
the phase transition a lot more are needed due to the 
critical slowing down. 

FigQ] is the phase diagram we obtained. Two tech- 
nical remarks are in place. First, to exam if the AFM 
to PM transition is continuous, we checked the inverse 
static spin susceptibilities at fc = (tt, tt, tt) which becomes 
very critical at the corresponding transition values of V 
[2^ . Second, the crossover temperature Tq is obtained 
using the saturated static homogeneous spin susceptibil- 
ity = 0,iO) C/Tq at low temperatures. [We used 
Curie constant C = 1/2 which is obtained in the high 
temperature limit.] An example of the saturation behav- 
ior is shown in the inset of FigQ] Note that To being small 
means the FM spin fluctuation becomes very strong. 
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FIG. 2: The local, XAA{iO), and nearest neighbor, XAB{iO), 
spin susceptibilities at /3 = 120 v.s. V. XAA(iO) is alway 
positive and shows a peak in the crossover region. XAB{iO) 
changes from negative in the AFM phase to positive in the 
Kondo phase. We multipled xab by the coordination number 
which reflects its contribution to the lattice spin susceptibility. 

To study the critical magnetic fluctuation around the 
QPT, we plot in Fig|2xAA(*0) and xab(^O). From the re- 
sult we see that XAB{iO) changes sign when V increases. 
This sign change is actually a special manifestation of 
a more general evolution of the RKKY correlation from 
being AFM to FM as hybridization is increased, due to 
the conduction electrons mediating the RKKY change 
from those near the Fermi surface to around the band 
bottom. Meanwhile, XAA{iO) is always positive and be- 
comes very strong in the crossover regime. This scenario 
should extend to T = and there would be a point at 
which XAsiiO) = and XAA(iO) becomes critical. This 
point would be the heavy Fermion QCP. Note that the 



static spin response being local does not mean a local 
electron self-energy. Actually, the non-local self-energy 
becomes stronger at lower temperatures p^. 
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FIG. 3: The static spin susceptibilities as functions of the 
temperature at (a) V — 0.26 and (b) V = 0.5. The fittings 
are given by x(T) = 7.399 ln(0.253/r) in (a) and x(T) = 
1.3141n(1.218/r) in (b). Note that according to Fig[T] V = 
0.26 is close to the QCP on the Kondo side. 

An important question is how the Kondo screening, 
which is local in space, becomes coherent in heavy 
Fermions. To this end, we study the evolution of the spa- 
tial correlation in the spin responses in Fig. |3| It shows 
that spin fluctuations are quite local in space down to 
T ~ 0.1 for y = 0.26 and T ~ 0.2 for V = 0.50. At lower 
temperatures, the FM spin susceptibility becomes domi- 
nant. This is similar to that observed experimentally in 
YbRh2Si2 |23||. Two remarks are in place. (1) the non- 
locality in X develops at lower temperature for V — 0.26 
than that for V = 0.50. This reflects the local nature of 
the spin fluctuation in the quantum critical region which 
extends to lower temperature as the QCP is approached 
and is consistent with Fig|2] (2) The logarithmic tem- 
perature dependence of x as shown in Fig. 13 is similar to 
those observed in many heavy Fermion compounds |[l|. 
A logarithmic temperature dependence is also found in 
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the total energy shown in Fig. 2J 
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FIG. 4: The total energy at V = 0.26. We fitted E{T) = 
-1.025 + 1.700r^ln(1.102/r). The corresponding specific 
heat coefiicient is 7(r) = 3.400 ln(0.6682/r). 

To conclude, using a two impurity DMFT, we studied 
the periodic Anderson model on cubic lattice at finite 
temperatures. We obtained the phase diagram which 
is consistent with the picture that both the Neel and 
Kondo temperatures vanish at the QCP. As the QCP 
was approached from the Kondo side, we found strong 
ferromagnetic spin fluctuations. From the sign change of 
the static nearest-neighbor spin susceptibility, we conjec- 
tured that the critical static spin fluctuation become local 



at the heavy Fermion QCP. We explored the crossover re- 
gion and observed logarithmic temperature dependences 
in the specific heat coefficient and spin susceptibility. 

Our results presented in this letter implies that a two 
impurity Anderson model combined with DMFT might 
serve as a minimal model in describing most of the ther- 
modynamics of the heavy Fermions. However, we should 
also keep in mind that there could exist other phases, 
like the superconducting phase , whose existence may 
require mechanisms which are more spatially extended 
than those describable in a two impurity model. 

Since the 4f or 5f orbital contributes to the local mo- 
ment physics in most heavy Fermion compounds, there 
are inevitablely many physical properties over simplified 
by the periodic Anderson model, like the orbital degener- 
acy and the subsequent crystalline field splitting. These 
are further complicated by the spin-orbital coupling, lat- 
tice frustration, disorder, hybridization with other bands, 
etc. Actually all these contribute to a much richer physics 
observed in experiments P| than what we have obtained. 
To this end, our current work can be considered as a use- 
ful guide to distinguish the "universal" heavy Fermion 
features from those specific to individual materials. 
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